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We make extensive simulations over a spin chain model that combines the frustrated J 1 -J 2 spin 
chain and the long-range nonfrustrated (—l)( r ~ 1 V _Q: decay interactions through the variational ma¬ 
trix product state method for both finite and infinite lengths. We study both the ground state entan¬ 
glement and phase diagram. We find that it is most entangled in the rotation invariant long-range 
ordered antiferromagnetic phase, where the entanglement scales approximately logarithmically. We 
determine the development of the Majudar-Ghosh point to a disorder line from entanglement. And 
we determine approximately the transition from the dimerized and incommensurate phase of the 
J 1 -J 2 model to a decoupled phase by studying spin correlation and the dimerization order parameter. 
Some implications for entanglement in systems with long-range interactions are stated. 


I. INTRODUCTION 

Quantum spin chains are the fertile ground to study 
strongly correlated quantum many body effects. One 
of the most studied models is the quantum Heisenberg 
model. For antiferromagnetic coupling (Ji > 0) the 
spin correlation of its ground state decays as 1/r up to 
a logarithm correction^ displaying quasi-long range or¬ 
der (QLRO). The model with a next-nearest neighbour 
added is known as the J 1 -J 2 or zigzag spin chain. It is 
frustrated when J 2 > 0. A dimerization transition occurs 
at J 2 /J 1 — 0.2411 after which it has a valence-bond-solid 
(VBS) orders and incommensurate spiral spin correla¬ 
tion emerges after the Majudar-Ghosh (MG) point at 

J 2 /J 1 = 0.5*4 

Beyond the next-nearest terms, the system can be built 
up with even long-range interactions (LRI). The models 
with power law decay of LRI coupling J r ~ (—l) r_1 r _a 
have attracted much attention^, from which, intriguingly 
true long-range ordered antiferromagnet (AFM) can be 
formed for small enough a£, even though it have been 
strictly ruled out from ID short-range rotation invariant 
models at even zero temperature. 

Recently, Sandvik proposed the combination of the 
J 1 -J 2 model and the long-range nonfrustrated terms for 
studying the interplay between then4 The Hamiltonian 
is 

N N 

i=i j=»+ 1 UJ 

J 2 = g , Jr ^2 = {-l) r ~ 1 r~°‘/n(N,a), 

where the normalization factor n(N, a) = 1 + r _a 

signifies that the sum of all the interactions (excluding 
r = 2) on one site add to unity, and it also ensures finite 
energy per site for infinite N when a < 1. The ground 
state phase diagram has been investigated by Sandvik^ 
using the exact diagonalization (ED) method, and later 
by Kumar and Soos^ using ED and other auxiliary meth¬ 
ods. Several phase boundaries have been accurately 
determined^, but it is still controversial in the regime 
with moderate frustration^^, where it should be related 


to the spiral state in the ab initio study of realistic metal¬ 
lic chains^. Besides, the ground state entanglement 
has not been considered in both works. In this paper, 
we restudy the model |T] using the matrix product state 
(MPS) approach for the ground state entanglement and 
phase diagram. As is established in recent years, entan¬ 
glement enriches the characterization of quantum phases 
and phase transitions^. Particularly in this model, it 
displays sudden drop along a first order phase transition. 
And a line segment of minimum entanglement marks the 
development of the Majudar-Ghosh point. This method 
also entails evaluation of spin correlation accurately for 
several hundred sites, so that we can directly demonstrate 
the incommensurate behaviour in the regime of moder¬ 
ate frustration. And we improved the phase boundary in 
this regime, which should resolve the controversy among 
the previous works. 

Another motivation of our work concerns the scaling of 
the entanglement entropy with subsystem size in quan¬ 
tum many body systems. In the last decade, much effort 
has been devoted to this problem^. Especially an area 
lawH of entanglement is conjectured, which states that 
for a local and gapped Hamiltonian the ground state en¬ 
tanglement entropy of a subsystem scales as the bound¬ 
ary area. This is remarkable, as it means that physical 
ground states are “slightly” entangled. In particular for 
ID systems it have been proved rigorously for gapped 
Hamiltonianii, and is shown to be violated mildly by 
logarithmic divergence for gapless and conformal invari¬ 
ant systems^. 

The area law is most generally attributed to the inter¬ 
actions being local. The entanglement structure in LRI 
systems is less known4 Especially, one wonders that 
should then the area law be severely violated? Some pro¬ 
gresses have been made in spin chains with LRI. In [13, 
it was shown that for the Lipkin-Meshkov-Glick model, 
which resembles the XY model but with infinite range in¬ 
teractions, the entanglement entropy scales at most log¬ 
arithmically. Lately, Koffel et. al.— studied a transverse 
Ising model with power law decay LRI, and showed inter¬ 
estingly that a gapped phase can even have logarithmic 
scaling of entanglement. That model is polarized, if the 
system is rotation invariant, stronger quantum fluctua- 
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tion and entanglement will present, which to our knowl¬ 
edge has not been considered so far. The model [T] being 
rotation invariant is ideal for addressing this question. 
Besides, with a tunable frustrating term, it facilitates to 
examine when LRI takes effect in increasing of entangle¬ 
ment. We find that it is most entangled in the AFM 
long-range ordered phase. But, remarkably, the scaling 
of entanglement can be still fitted approximately with 
logarithm functions. This indicates that, in contrary to 
one might expected, the area law should be not severely 
violated in this system. 

This work is organized as follows. In Sec [II] we briefly 
introduce the numerical methods. See lllll presents the 
simulation results. In Sec. IIII Al we show the ground 
state phase diagram, half chain entanglement entropy, 
and the development of the Majudar-Ghosh point. Some 
numerical pitfalls and partial solution to them are also 
manifested. In Sec lIII Bl we discuss the scaling of en¬ 
tanglement with subsystem sizes. In Sec lIII Cl we study 
the phase boundary at moderate frustration. Finally, we 
conclude in Sec lIVI with some implications stated. 


II. METHOD 

It is conventionally difficult to simulate systems 
with LRI using density matrix renormalization group 
(DMRG)^ or MPS alike methods. Crosswhite et. al3£ 
have made a major progress for LRI systems with power 
law type decay by approximately encoding the Hamil¬ 
tonians with matrix product operators (MPO), so that 
the computation cost is reduced considerably. We use 
this long range MPO method to represent EqJT] for fi¬ 
nite and infinite systems^ in the parameter region of 
(a, g) 6 [0.7, oo) x [0,1.0] (we refer to [2^ and I 22 I for de¬ 
tails). A minor adaption for the J 2 term is needed. We 
write Eq[l] in an equivalent form, 



FIG. 1. Ground state phase diagram for EqQ] At a = 00 , the 
model degenerates to the J 1 -J 2 chain, whose phases undergoes 
phase transitions under LRI with decreasing a: the QLRO 
phase to a long range order AFM phased; the commensurate 
VBS phase terminating at a multi-critical point at around 
a = 1.7 and g = 0.41 and continued by a line of first order 
transition (thick solid line)^; the incommensurate VBS phase 
to a phase decoupled for odd and even sublattice. 

DMRG by White) is used to study infinite systems. This 
method exploits the translation invariance, such that 
the computation effort is reduced and boundary effect 
is avoided. It is implemented with a 4-site unit cell, from 
which an infinite MPS (IMPS) representation can be re¬ 
constructed after convergence^ for measuring physical 
observables. The maximal D used is 1000 for generat¬ 
ing a well converged fixed point with truncation error 
restricted to smaller than le-6, while at some point we 
also use iDMRG to generate a finite open chain with even 
larger D. 

For a ground state wave function | ip) on a finite open 
chain [1, N] of length N, divided into a subsystem [1, L] 
and its complement [L + 1, N], we measure the entangle¬ 
ment entropy 

S =-Tr(p L \n(p L )), (3) 


H = H a + H b 

**-!ob*E i (-l) li -‘ | - 1 |j-r“SS; ,,, 

i— 1 j=i +1 \*J 

N -2 

^ = (3+2^1) E Wm. 

2=1 

Now H a has uniformly power law decay interaction (with 
alternation signs), so that the long range MPO method 
is applicable. The short-ranged H b can also be encoded 
in a MPO easily. Then the entire Hamiltonian is the sum 
of the two MPOs22. 

We use the variational MPS (VMPS) algorith m 22 ' 23 to 
simulate the ground states for finite open chains with N 
ranging from 16 to 100 and truncation dimension D up 
to 520. It is implemented with the 1-site algorithm and 
density matrix correction? 4 that reduces the chance be¬ 
ing stuck. The quality of the variational ground state is 
gauged by the average variance v = (( H 2 ) — ( H) 2 )/N , 
kept smaller than le-4 for the hardest case. The iDMRG 
algorithm?^ (not exactly the conventional infinite size 


with pl = Tr^ +1 jv](|'0)(' ! /’l) the reduced density matrix 
for the subsystem. We also measure other quantities that 
will be defined later. 


III. SIMULATION RESULTS 
A. ground state entanglement and phase diagram 

We plot a schematic ground state phase diagram in 
Fig. [T]and show the distribution of entanglement on the 
parameter plane in Fig[5] The entanglement is generally 
higher as a reduces; rather high in the top left corner 
while low in the bottom center. These already give a 
rough profile of several of the phase boundaries. In the 
center of the phase diagram is a A-shaped phase bound¬ 
aries, with the dimerized and commensurate phase be¬ 
neath them. Sandvib^ has successfully determined the 
development of dimerization point of the J 1 -J 2 model to 
a multi-critical point at around (1.7, 0.41), and then to 
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FIG. 2. Distribution of the half chain entanglement entropy 
Sn /2 on the parameter domain (a, g) £ [0.7,3.0] x [0,1.0] 
divided into 24 x 50 points for system size N = 60. 



FIG. 3. Dependence of half chain entanglement Sjv /2 on g 
for several values of a with N = 100. Inset displays in great 
detail of the lines around g = 0.4 for a = 1.6, 1.7 and 1.8. 

a first order phase transition, using level crossing. Be¬ 
low, we study these phase boundaries again but from an 
entanglement perspective, and we determine the devel¬ 
opment of the Majudar-Ghosh points as well. 

FigO shows in particular dependence of entanglement 
on g for several values of a. At a = oo, there are two 
turning points for the curve: The first one is related to 
the dimerization transition point at g c — 0.2411 (where a 
gap opens and entanglement drops), but it is difficult to 
locate g c accurately from entanglement; The other one 
is the minimum of entanglement just at the Majudar- 
Ghosh point gMG = 0-5- The two points are smoothly 
connected and approaches as a reduces until the point 
(1.7,0.41) (see inset of Fig[3j), indicating shrinking of 
the VBS(comm.) phase and finally terminating at that 
point. That multi-critical point is in agreement with the 
ED result obtained by extrapolation^. The trajectory 
of gMG as a changes appears as a valley in the entan¬ 
glement plane of Fig(21 where it should also have mini¬ 
mum correlation length, and thus can be thought of as 
a disorder-line^, separating phases with commensurate 
and incommensurate correlation on either side. After the 
multi-critical point, entanglement becomes discontinuous 
displaying a sudden deep drop, which clearly marks the 



FIG. 4. (a) Dependences of variance (open square) and mid¬ 
dle chain entanglement (open circle) on a shown together for 
comparison for fixed g = 0.0 and N = 100, D = 520. (b) 
Ground state energy as a function of g at a = 1.0 obtained 
from different algorithms. Open triangle: VMPS initiated 
from iDMRG or random state; open square: VMPS with 
proper initial state; cross: iDMRG. Inset shows extrapolation 
of position of peak of energy using N ranging from 16 to 100. 
The fitting line is g = 1.0 /N + 0.339. (c) Middle chain en¬ 
tanglement as a function of iDMRG steps at a = 1.0, g = 0.0 
with D = 1000. 


first order phase transition from the AFM phase to the 
decoupled phased. 

The rather high amount of entanglement in the AFM 
phase poses considerable difficulties to the MPS method, 
since the computational effort for it scales exponentially 
with entanglement^. Here we would like to elaborate 
on these difficulties. Fig0][a) shows dependence of both 
average variance and entanglement on a at g = 0.0. This 
serves as a benchmark of the accuracy viable. One can 
see that variance increases radically with decreasing a. 
These restrict us to cr > 0.7 for the variance smaller 
than le-4 for the maximal length N = 100 and largest 
D = 520 used. And we find much more sweeps (around 
10 times) needed for convergence for small a. Further 
more, there is metastable state issue to the left of the first 
order transition point. As shown in FigUKb), it is prone 
to get stuck on an excitation level which should have less 
entanglement than the ground state, if one uses random 
state or iDMRG for an initial state. This leads to a wrong 
position of the peak of energy (the transition point) com¬ 
pared with ED^. For a given N, larger D can shrink the 
region being stuck, but soon become unpractical. A two- 
site algorithm with density matrix correction won’t solve 
it either. It turns out a nice solution is to provide a better 
initial state, e.g., use the state of a smaller g as the input 
of larger g close to the right of the boundary (see j2Sj for 
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alternative ways such as adding a pinning term for fixing 
this). In this way, the peaks for each lengths are unam¬ 
biguously determined and the extrapolated value of the 
transition point is g c = 0.339 (see inset of Fig@fb)). As 
for the infinite algorithm, the metastable issue is more 
severe. It is stuck in a wider range, which, we however 
haven’t found a way to avoid. For g close to 0 the energy 
still deviates with VMPS. This is not because of getting 
stuck but is a convergence problem due to too fast grow¬ 
ing of entanglement and at the same time relatively slow 
convergence of energy. As shown in Fig0|c), entangle¬ 
ment suddenly drops after around 30 iteration steps (120 
chain length) if D = 1000 is kept not increased and even¬ 
tually converged to a wrong fixed point. One could stop 
iteration before the drop (this is where the data points of 
energy we adopted), but energy and other quantity are 
far from convergence. In all we find good convergence of 
VMPS for the parameter range studied, while iDMRG 
has either metastable state or convergence problems in 
the AFM phase. 


B. scaling of entanglement under LRI 

As shown above, the highest amount of entanglement 
is found when frustration is zero and a is small. At 
g = 0.0, according to previous ED^ and quantum Monte 
Carlo^ studies, the system undergoes a continuous phase 
transition from the QLRO to AFM phase at a c ~ 2.22. 
Recall that for a conformal invariant system with open 
boundary condition, the entanglement scales logarithmi¬ 
cally with subsystem length— 

5~|ln(L), (4) 

where c in the prefactor is identified with the central 
charge of the relevant field theory. The QLRO phase be¬ 
longs to this category and has c = 1. We expect changing 
of behaviour of S(L) around a c , and wonder to what ex¬ 
tent the entanglement in the AFM phase with strong LRI 
violates the area law. 

To this end, we focus primarily on the line of g = 0.0 
and use both VMPS and iDMRG to simulate various 
sized systems. The former, being free of environment 
error— and better controlled, is used for N < 100, while 
the later having reduced computation cost, is used for 
100 < N < 200 and validated by the consistency with 
small chains. We first present S(L) for the largest size 
N = 200 for each a in Fig|5la). One can see remarkably 
that, for all a values S scales approximately logarithmi¬ 
cally with L (for L not too close to the chain center), 
and that S(L) increases clearly faster for smaller a. To 
quantify and compare them, it is tempting to use Eq0] 
to fit each curve and extract a value of c. We may call 
c generally an “effective central charge”, as the AFM 
phase is not conformal invariant^. As shown in the inset 
of Fig(5ja), c is near to 1.0 for a < 2.5, after that, it in¬ 
creases quickly with a and reaches a large value. This be- 




FIG. 5. (a) Bipartite entanglement entropy S as a function of 
subsystem length L for each a with g = 0.0. Only bipartition 
on even bonds (even L) are drawn. Solid lines are fitting 
to S — | ln(L) + const for 1 < L < 75. Inset shows the 
extracted c vs. a. (b) Three groups of curves at (1.0, 0.0) 
(open circle), (1.0, 0.7) (open square) and (3.0, 0.0) (open 
triangle) each shows dependence of S on L for different system 
sizes N ranging from 32 to 200. L is restricted to no larger 
than | N, i.e. not close the chain center. Solid lines are fitting 
to the function as in subgraph (a). Inset shows c extracted 
from fitting of each group of curves as a function of N. 


haviour is in overall agreement with the transition point 

G? c . 

The above extracted values of c in the AFM phase are, 
however, not quite validated. For one thing, closely ex¬ 
amining the fitting, one finds small deviations from per¬ 
fect logarithm: S(L) seemingly increases slightly faster 
for larger L. For another, and more significantly, S(L) 
actually has a salient dependence on the total system 
size, as a result the c values are only specific to N = 200. 
To manifest this finite size effect on S(L), we compare 
three points (1.0, 0.0), (1.0, 0.7) and (3.0, 0.0) which 
are respectively representatives of the AFM, decoupled, 
and QLRO phases. For each point we plot S(L) for var¬ 
ious system sizes, as shown in Fi^SJb). Since they all 
have (approximately) logarithm divergence, the values of 
c are extracted for each N, and shown in the inset of the 
graph. One can find that a clear dependence of S(L) and 
c on N is unique to the AFM phase. (In the decoupled 
phase, c are close to 2.0, which is expected. Because, as 
will be shown later, it is a system of two approximate 
Heisenberg chain, whose central charge is just the sum of 
each one’s.) The dependence of S(L) on N implies that 
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the maximal size N = 200 should be not enough. We will 
try to give an extrapolated result for it later, but below 
we first try to interpret the results obtained. 

So far in the above we observed that for all a values, 
S(L) have approximately logarithm dependence, while 
the slope dS/d{ln(L)) i) apparently increases for smaller 
a, ii) slightly increases for larger L and iii) increases 
for larger A. We give a naive explanation for these 
from a valence-bond description of the singlet ground 
state of antiferromagnet. A valence-bond is a singlet 
pair - 7=1 f iij ~ where i and j are on opposite 

sublattices of of a bipartite lattice. It is known that, a 
SU(2) singlet ground state can always be represented in 
a valence-bond basis which is all possible covering of sin¬ 
glet pairs on the chain [1, AT] for even TV. Each singlet 
pair is maximally entangled, and the value is ln(2). Thus 
bipartite entanglement can be measured as the number 
of bonds cut by the bipartition times ln(2)21. This give a 
appealing geometrical interpretation of the ground state 
entanglement. For the Heisenberg model, Si ■ Si +1 favors 
forming of singlet between spin i and i + 1 , but many 
body effects eventually leads to forming of complex dis¬ 
tribution of configurations of valence-bond^I - including 
pairs apart in arbitrarily long distances. Note it is the 
long distance entangled pairs that lead to divergence of 
entanglement with increasing L. Notably, for the un¬ 
frustrated regime of the Hamiltonian [T| long-range terms 
Jj-iSi ■ Sj (we may call them bonds of the Hamiltonian) 
favor directly long distance singlet pairs. It is then rea¬ 
sonable to assume that the number of valence-bond in the 
ground state across [1, L] and [L + 1, TV] is in a way posi¬ 
tively correlated with the sum of the strength of all bonds 
of the Hamiltonian across them, i.e. with the quantity 

J = y (j - i)~ a . ( 5 ) 

(Note that we can safely ignore the normalization in the 
Hamiltonian when g = 0.) Based on this, we may under¬ 
stand the properties just stated. For i), J increases with 
smaller a, so entanglement increases faster as well. (But 
this does not explain why they are close to logarithm for 
all the a); For ii) and iii), J sums over L x (TV — L) 
(approximately L x TV for large TV and small L ) bonds, 
while in contrast, for the QLRO phase which is essentially 
short-ranged, there is always only one bond connects the 
two regions irrespective of the position of the cut or the 
system size. This helps to understand why in the AFM 
phase the approximate logarithm function S(L ) increases 
even faster when L or TV expands, while it does not in 
the QLRO phase. 

To get a glance of S(L) in he thermodynamical limit, 
we try to extrapolate the values S(a, L, TV) in TV for each 
fixed a and L. Hereafter we make the dependence of S 
on a, L and TV explicit, leaving fixed g = 0.0 implicit. 
In FigH>ta) we fix L = 12 and show for different a the 
dependence of S(a, L , TV) on 1/TV. In FigOb), we fix a = 
1.0 and show for different L the dependence of S(a, L , TV) 
on 1/TV. They show that the extrapolation is plausible. 





FIG. 6. (a) S(a,L,N) versus inverse chain length for dif¬ 
ferent a and fixed L = 12. Solid lines are fitting to 
S(a,L,N) = — k/N + S(a,L, 00 ) which is a linear function 
of 1/A. (b) S(a,L,N) versus inverse chain length for dif¬ 
ferent L and a = 1.0. Solid lines are fitting to the same 
function as in subgraph (a). The system sizes A used are re¬ 
stricted to [Nmin, 200], with Nmin depending on L, such that 
§A m in > L, namely, L should not close to the chain center, 
(c) Extrapolated S(a, L, A) value to A = 00 as a function of 
L (L < 64) for several a, obtained by curve fitting as above. 
The lines are guide for eyes. 

For different L and a there is always a 1/A relation 

S(a,L, N) = +S(a,L,oo), ( 6 ) 

but with a non-travail coefficient k depending on a and 
L. Two profiles of k(a, L) for fixed L = 12 or a = 1.0 are 
shown in the insets of Fig]6]]a) and (b). The finite size 
effect is more salient for smaller a or larger L. This is 
in agreement with the above argument from the valence- 
bond description. 

Then extrapolated S(a,L,oo) for more a and L com¬ 
binations are shown in FigH^c). Note that it should 
be understood as a result for semi-infinite chain [l,oo). 
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FIG. 7. Absolute value of the intra sublattice correlation C(r) 
for each a at g = 0.7. C(r) are obtained by evaluating Eq[7| 
upon IMPS from iDMRG simulations with D up to 960. Inset 
shows C(r) at a = 1.4 and g = 0.7 multiplied by y/re r A with 
£ = 44.7 using iDMRG and D = 640. 

One can see that each curve is still close to a loga¬ 
rithm function, with the slop dS/d(ln(L)) slightly in¬ 
creases with larger L. The maximal L evaluated is 
rather restricted, it is not very clear whether dS/d(ln(L)) 
keeps increasing slowly or approaches a constant for even 
larger L. For the former, the function form can be e.g. 
S = a\ ln(L ) 2 + 02 ln(L) + 03 , while for the later, it can 
be e.g. S = ai ln(L) + 0 , 2 /L + 03 . We have verified that 
both can give a good fit for the curves (not shown), other 
than a simple logarithm function. But in view of the 1 /N 
correction in Eq| 6 l and that in Figjsjb) curves for differ¬ 
ent N are all close to logarithm functions, we are prone 
to the later. 


C. phase boundary at moderate frustration 

We next turn to the right part of the phase diagram. 
Focusing on one line g = 0.7, we measure the spin corre¬ 
lation 

C(r) = (S*S? +r ) (7) 

and see how it changes as a reduces. Here the spin chain 
is considered to be divided into odd and even sublattices, 
as is usually did for the J 1 -J 2 model. Fig(3 shows corre¬ 
lation between spins within a same sublattice (even r). 
One can see that for large a values C(r) decreases ex¬ 
ponentially with jumps in the curve. The jumps signify 
the incommensurate behaviour: Following the treatment 
of White and Affleck^ of the J 1 -J 2 model, we multiply¬ 
ing C(r) e.g. at a = 1.4 and g = 0.7 by y / re 7 '^, then 
the sinusoidal modulation is clearly seen in the inset of 
the graph, where the correlation length £ = 44.7 is cho¬ 
sen such that the beats of the amplitude are as flat as 
possible. By evaluating C(r ) for r up to 1000, we only 
find jumps for a > 1.2. For a < 1.2, C(r) displays alge¬ 
braical decay (for reference C(r) ~ 1/r 7 with 7 = 1.18 




FIG. 8. Left and right shows respectively dimerization d and 
minus of inter sublattice correlation for r = 1 , measured upon 
IMPS, as a function of iDMRG truncation dimension D for 
several (a, g ) tuples. 


at a = 1.0). This indicates a possible critical value a c of 
a continuous phase transition at roughly 1 . 2 . 

In previous studies with short chains, Sandvik^ esti¬ 
mates that the transition should be from the VBS phase 
to a coexisting VBS+QLRO( 7 r/ 2 ) phase, while Kumar 
and Soosi£ predict that the transition is from VBS phase 
to a decoupled phase. We agree with Kumar and Soos 
that for small a there should be no VBS order and 
is essentially decoupled, but from an independent and 
more direct way: In fact the correlation shown above 
already implies no VBS order for a < a c , because the 
presence of incommensurate behaviour is related to a 
finite dimerization^. We also show directly the dimer¬ 
ization order parameter^ d = {SfS* +1 ) — {SfS- : _i), 
for several (a, g) tuples in Figj 8 ](a). One can see that, 
for a = 2.0 > a c , d are clearly nonzero; While for 
a = 1.0 < a c , d are very small and appear to vanish 
in the limit D —> 00 . In addition, the inter-chain correla¬ 
tion is very small for a < a Cl as shown in Figj 8 ])b). (Here 
it suffices to consider C(r) for r = 1, since C{r ) further 
decays for larger odd r .) 

The transition points for other g can be determined 
likewise and we find that a c increases with g , and it 
should be that a c —> 00 as g —> 00 . This gives a approx¬ 
imate phase boundary between VBS(incomm.) and the 
decoupled phase in Fig(T] Note crucially that the starting 
point of the phase boundary is distinct from both Kumar 
and Soos’s or Sandvik’s result, in which, it starts at a 
unique multi-critical at around (1.7, 0.41). The signifi¬ 
cance of this is that there indeed can be direct transition 
from AFM long-range order to VBS order—. 

The phase transition can be understood as follows: 
At moderate or large < 7 , as a reduces, the couplings J r 
for r ^ 2 (including J{) all become very small due to 
a large normalization factor, and eventually the next- 
nearest neighbour term gSi ■ S )+2 dominates. The later 
induce a background of antiferromagnetic order on ei¬ 
ther sublattices, this amounting to a period 4 structure 
in the entire lattice. It is not difficult to see that, the rest 
enormous but small couplings J r ^2 now have conflicting 
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signs against this background period, and their effects 
should be largely smeared out. So for mall enough a it 
is decoupled into nearest-neighboured Heisenberg model 
with coupling g and with long-range terms as perturba¬ 
tions. Note that for the J\-Ji model (a = oo), White 
and Affleck^ used field theory to predict that there is ex¬ 
ponential small gap and dimerization for arbitrary large 
g , except for g = oo where it decouples exactly into two 
Heisenberg chain, and supported it by DMRG. While 
our arguments above essentially states that, at small a, 
the spin chain can be decoupled for modest g. Never¬ 
theless, the evaluation of C(r) is still not very long and 
the decay of d with D is somewhat slow, so we are still 
not completely sure whether there are incommensurate 
modulation with very long period or dimerization (and 
also spin gap) should be exponentially small but nonzero 
even as a approaches 0 , which is very difficult to confirm 
numerically. A field theory study for small a may be 
desirable as well as that for the J 1 -J 2 model at large g. 


IV. CONCLUSION 

In summary, we studied the frustrated spin chain with 
long-range interactions using the matrix product state 
approaches. We found that it is most entangled in the 
rotation invariant AFM long-range ordered phase, where 
the entanglement scales approximately in logarithmic 
form. But we miss a complete understanding for the 
logarithmic scaling. The maximal systems size studied 
is not large and a is still not close to 0. It is worth¬ 
while to check the persistence of the logarithmic form 
in the asymptotic scaling of entanglement for large sub¬ 
system size and smaller a. in future works. We also stud¬ 


ied correlation and dimerization for moderate frustration 
and determined an approximate boundary for the tran¬ 
sition from the dimerized and incommensurate phase of 
the J 1 -J 2 model to a decoupled phase. 

Our work implies that in an antiferromagnetic spin 
chain with LRI, frustration terms usually hamper long- 
range entanglement. Consider a series of power law de¬ 
cayed models H({s}) = r s(r)r~ a Si ■ Si +r , with the 
signs s(l) = +1 while s(r > 1 ) £ {+ 1 ,- 1 } indetermi¬ 
nate. At small a, their entanglement should be bounded 
by the unique nonfrustrated one i.e. s(r) = (—l) r ’ _1 , 
where perfect AFM long-range order is formed. 

The scaling of entanglement shown in this work and 
previous work s 17 ’ 18 indicates that, although LRI can usu¬ 
ally increase entanglement, but it not necessarily leads to 
large entanglement and severe violation of area law. The 
preservation or (severe) violation of area law may be a 
joint effects of interaction range, symmetry constrains 
(e.g. translation invariance), ground state degeneracy, 
and etc, which remains to be clarified. 
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